Draft version February 26, 2013 

Preprint typeset using I^TgX style emulateapj v. 5/2/1 1 



(N 

o 



m 



O 
U 

O 



OO 

rn 

O 



% 



PROBING THE INFLATON: SMALL-SCALE POWER SPECTRUM CONSTRAINTS 
FROM MEASUREMENTS OF THE CMB ENERGY SPECTRUM 

Jens Chluba', Adrienne L. Erickcek'-^ and Ido Ben-Dayan'-^ 
Draft version Februarys 26, 2013 

ABSTRACT 

In the early Universe, energy stored in small-scale density perturbations is quickly dissipated by Silk- 
damping, a process that inevitably generates ju- and y-type spectral distortions of the cosmic microwave back- 
ground (CMB). These spectral distortions depend on the shape and amplitude of the primordial power spectrum 
at wavenumbers k < 10"*Mpc"'. Here we study constraints on the primordial power spectrum derived from 
COBE/FIRAS and forecasted for PIXIE. We show that measurements of // and y impose strong bounds on 
the integrated small-scale power, and we demonstrate how to compute these constraints using ^-space window 
functions that account for the effects of thermalization and dissipation physics. We show that COBE/FIRAS 
places a robust upper limit on the amplitude of the small-scale power spectrum. This limit is about three orders 
of magnitude stronger than the one derived from primordial black holes in the same scale range. Furthermore, 
this limit could be improved by another three orders of magnitude with PIXIE, potentially opening up a new 
window to early Universe physics. To illustrate the power of these constraints, we consider several generic 
models for the small-scale power spectrum predicted by different inflation scenarios, including running-mass 
inflation models and inflation scenarios with episodes of particle production. PIXIE could place very tight 
constraints on these scenarios, potentially even ruling out running-mass inflation models if no distortion is de- 
tected. We also show that inflation models with sub-Planckian field excursion that generate detectable tensor 
perturbations should simultaneously produce a large CMB spectral distortion, a link that could potentially be 
estabUshed by PIXIE. 
Subject headings: cosmic microwave background - theory - observations - inflation 



1. INTRODUCTION 



Co s mological in flation dAlbrecht & SteinhardtllI982l: iGuthI 
II98II: lLindelfT982b provides a commonly accepted explana- 
tion for both the Universe's homogeneity and the origin of the 
initial curvature perturbations that seeded the growth of struc- 
ture. Inflation cannot be considered a complete theory, how- 
ever, until we understand the inflaton: the field that drove an 
epoch of accelerated expansion in the early Universe. Fortu- 
nately, the statistical properties of the initial density perturba- 
tions offer a wealth of information about inflationary physics. 
In single field inflation, we can in principle reconstruct the 
inflaton potential if the primo rdial power spectrum is known 
at all scales (e.g. iLidsev et al.iil997i) . However, the limited 
range of scales probed by the CMB and large scale structure 
(LSS) does not provide sufficient information to discriminate 
between many inflation models. Finding additional ways to 
measure the primordial power spectrum outside this range of 
scales will greatly enhance our ability to constrain the infla- 
ton's potential and its trajectory during inflation. In this work, 
we investigate how spectral distortions of the CMB caused 
by the dissipation of energy stored in small-scale density per- 
turbations can provide a new probe of inflation by extend- 
ing our knowledge of the primordial power spectrum from 
k ^ 1 Mpc"' to about lO'^Mpc"'. 

The simplest models of inflation predict a power spec- 
trum parameterized by a nearly constant, slightly red spec- 
tral index. More complicated inflationary models can 
leave distinctive imprints in the primordial power spectrum. 



Multi-field infl ation can produce primordial power spec- 
tra with steps (Silk & Tu rned II9871 [Polarski & Starobinsk^ 
1992; Adams et al. 1997|) or oscill a tions ( Achilc arr o et al. 
201 U iKobavashi & Takahashil 1201 It ICespedes et al] |20I2[k 






These features in the primordial power spectrum may also 
be generated during single-field inflation by d iscontinuities, 
kinks, and bumps in the infl aton pote ntial (Salocek et al. 
J^89; Starobinskii 1992; Iv anoveta.| _ 1994; Starobinsk 
lT998t IHunt & Sarkari 120071 IJoyetal.l 12008). The primor- 
dial power spectrum may also contain information about how 
the inflaton interacts with other fields; for instance, parti- 
cle production durin g inflation leaves a bump in the primor- 
dial power spe ctrum Ichung et al.ll2000l: iBarnaby et"al1l2009l 
lBarnabyl l20I0l) . Finally, several inflationary models predict 
enhancement of the small-scale perturba tions that are gener- 



ated during the l ater stages of inflation (Randall et al. 1996; 
Steward Il997bt iCopeland etd] fT998t iCovi & Lvth .1999 



2001 
L^ 



Covi et alJII999l;lMartin et al.ll2000tlMartin & Brandenbergeri 



Ben-Davan & Brustein' '2010t iGong & Sasakil l2oTlt 



2011a; Bugaev & Klimai 201 IJ). 



The CMB temperature fluctuations provide a precise mea- 
surement of the primordial power spectrum on large scales, 
corresponding to wavenumbers 10'^ Mpc' < A: < O.I Mpc ' 
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(iPearson et al.i20()3l; lReichardt et al.l20 09tlBrown e t al.l2009l; 
iLarson et al.ll201 iF lHlozek et al.l l20I lb . Luminous red galax- 
ies and galaxy clusters probe the matter power spectrum on 
similar scales (0.02 Mpc~' < k < 0.7Mpc~ '; [Reid et alJ 
20TO [ rV^khTinin et all 120091; [Tinker et al.l 120121; [Sehgal et ail 
201 II) . while the Lyman-a forest reaches slightly smaller 
scales (0.3 Mpc' < k < 3 Mpc' ; [McDonald et al.1 12006[). 
All these observations indicate that the primordial power 
spectru m is nearly scale-invariant with an amplitude close to 
2x10'' (iTegmark & Zaldarriagal2002llNicholson & Contaldil 
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2009 [lKomatsu et alJlIoTTtlDunklev et aI][20Tll :[Keis"ler et alj 



20Tlt iHlozeketalJ I20TTI: iBird et alJ l20Tlh . There is no 
evidence of fe atures in the prim ordial power spectrum 
on these scales ([Kinney etalJ 1200 8: Mortonson et al. 200^ 



Barnabv & Huana 2009; Hamann et al. 2010; Peiris & Verd 
2010; Dvorkin & Hu 2010; Bennett et al. 201 1; Benetti et al 
20n||bvorkin & Hu 201 ij 



Our knowledge of the primordial power spectrum on 
smaller scales is far more limited; we only have upper 
bounds on its amplitude for k > 3Mpc '. One of 
these upper bounds is derived from the limits on spec- 
tral distortions in t he CMB. I t was long understood that 
the Silk-damping (ISiUd Il968h of primordial small-scale 
perturbations causes energy r elease in th e early Universe 
jSunvaev & Zeldovich| \l92M iDalvl 119911; iBarrow & Colesi 
[199I; Hu et al. 199#^. This gives rise to small spectral 
distortions of the CMB spectrum that directly depend on 
the shape and amplitude of the primordial power spectrum. 
Modes with wavenumbers 50Mpc"' ^ ^ ^ lO^Mpc"' dis- 
sipate their energy during the ju-era (redshift 5 x 10^ < 
z < 2 X 10^), producing a non-vanis hing constant residual 
chemical potential at high frequencies dSunvae v & Zeldovich 
1970bt IZefbovich et alj [19721; llllarionov & Sunvaev .1974: 



Buriganaetal.lll991l:lHu&Silklll993ah . while modes with 
k < 50Mpc"' result in a y-distortion. The latter is also well 
known in connection with the SZ-efFect of clusters of galaxies 
(IZeldovich & Sunyaevlll969h . By accurately measuring the 
CMB spectrum one can therefore place robust upper limits on 
the possible power at small scales since the physics going into 
the production of these distortions is well understood. 

Very precise measuremen ts of the CMB spect r um were ob- 
tained with COBE/FIRAS (iMather et all 1 19941; iFixsen et alJ 
[1996), limiting possible deviations from a blac kbody to /i < 
9 X 1 0'^ and y < 1.5 x 10"^ at 95% confidence (iFixsen et al.l 
[l99a). At lower freque ncies, a similar limit on fi was recently 
obtained by ARCADE (ISeiffert et al.ll20 Ih. and ^ < 6 x 10'^ 
at V :^ 1 GHz is derived from TRIS (IZannoni et al.l l2008t 
iGervasi et al.ll2008h . For power spectra with constant spectral 
index, «s, and normalization fixed at CMB scales, the mea - 
surements of COBE/FIRAS imply «s < 1.6 (iHu et alj|1994h . 
but this limit is model-dependent. For instance, a small neg- 
ative running ^ «nin, o f the spectral index weakens this bound 
significantly (JKhatri e t al. 201 1; Chluba et al. 2012). 

Here we generalize the COBE/FIRAS limits on spectral dis- 
tortions by directly converting them into a bound on the total 
perturbation power at small scales. Depending on the par- 
ticular inflationary model, this translates into constraints on 
different model parameters; the conversion can be obtained 
on a case-by-case basis. Also, the re cently proposed CMB 
experiment PIXIE ( iKogut et alj|201l[) might be able to de- 
tect distortions that are ~ 10^ smaller than the upper limits 
given by COBE/FIRAS. At this level of sensitivity, PIXIE 
is already close to what is required to detect the distortions 
arising from the dissipation of acoustic modes for a power 
spectrum with «s = 0.96 and no running all the way from 
CMB-anisotropy s cales to k =^ 10"^ JChluba & Sunvaevll20l"2[; 
iKhatri et alll2Cil It IChluba et aT]|2012l) . Such an improvement 
could rule out inflationary models with additional power at 
small scales, as we discuss here in more detail. Conversely, 
any detection of spectral distortions implies that either the 
power spectrum is enhanced on small scales, contrary to the 
predictions of the simplest inflation models, or an alternative 
mechanism generated CMB spectral distortions in the early 



Universe (e.g. particle decays). 

The only other upper bounds on the amplitude of the small- 
scale primordial power spectrum are derived from the absence 
of primordial black holes (PBHs) and ultracompact miniha- 
los (UCMHs), which are dense dark matter halos that form at 
high redshift (z ^ 1000). Both PBHs and UCMHs form in 
regions with large primordial overdensities; an initial over- 
densi t y of 6p/p > 0.3 is require d to form a PBH (ICarrI 
119751; (Niemever & Jedam"zigil999[) . while UCMHs form in 
regions where 6p/p > lO""* when they enter th e Hub ble hori- 
zon (Ricotti & Gould 2009; Bringmann et al.l 1201 Ih . There 
are numerous con straints on the number density of PBHs; 
iJosan et all ( 120091) showed that these constraints imply that 
the amplitude of the primordial curvature power spectrum is 
less than 0.01-0.06 over an extremely wide range of scales 
(0.01 Mpc' < k < lO-^Mpc'). Even though PBHs pro- 
vide only a weak upper bound on the small-scale ampli- 
tude of the primordial power spectrur n, they have usefully 



constrained inflat i onary models (e.g.. ICarr & LidsevI 



Leach et al.1 [20001: iKohri etal.112008) ; 'Peiris & Easthei 



1993 



2008 



Josan & Greenl2010al;lLvthl201 lbl: lBugaev' & Klimaii201 Ibl) 



Since UCMHs form in lower density regions than those 
that produce PBHs, limits on their abundance can pro- 
yide tighter constraints on the primordial pow er spectrum 
( IJosan & GreenI l2010bl: iBringmann et al.l 120111) . Unfortu- 
nately, all current limits on the number density of UCMHs 
rely on the assumption that they emit gamma rays from 
the annihilation of dark matter par t icles within their high 



densit y centers (Scott & Sivertsson 2009; Berezinsk v et al. 
i2010j :iLacki & Beacom 2010; Yang et al. 201 la b ; Yang et al 



1201 U [Zhang 2011; Bringman n et al.ll201 ll). If dark matte r is 
a self-annihilating thermal relic. IBringmann et al.l (12011b re- 
cently showed that the Large Area Telescope on the Fermi 
Gamma-Ray Space Telescope (lAtwood et al.l2009l) places the 
strongest constraint on the UCMH abundance; this limit im- 
plies that the amplitude of the primordial curvature power 
spectrum is less than 2 x 10"^ to 2 x 10"^ for modes with 
lOMpc"' < k < lO^Mpc"'. If dark matter does not self- 
annihilate, then UCMHs can only be detected gravitationally. 
In this c ase. [ETet al. (2012,) re cently showed that the Gaia 
satellite dLindegren et al.|[2012l) will be able to detect astro- 
metric microlensing by UCMHs and that a null detection of 
UCMHs by Gaia would constrain the amplitude of the primor- 
dial power spectrum to be less than 10"^ for k ^ 3500 Mpc"' . 
CMB spectral distortions probe the amplitude of the pri- 
mordial power spectrum in a very different manner than PBHs 
and UCMHs. First of all, the physics underlying the com- 
putation of the CMB spectral distortions is very well un- 
derstood, while, for example, the constraints derived from 
UCMHs depend on unknown properties of the dark matter 
particle: its mass, the abundance of its antiparticle, and its an- 
nihilation cross section. Second, since they arise from over- 
dense regions, PBHs and UCMHs probe the high-density tail 
of the probability distribution function for density perturba- 
tions. The likelihood of forming a PBH or an UCMH is there- 
fore highly sensitive to deviations from Gaussianity that en- 
hance or suppress the abundance of high overdensities. In 
contrast, the CMB spectral distortion is determined by the 
total energy stored in density perturbations and is therefore 
less sensitive to the precise form of the probability distribu- 
tion function. Furthermore, the expected number density of 
PBHs or UCMHs of a particular mass is determined by the 
mass variance in the sphere with radius R that formed the PBH 
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or UCMH, which is computed by convolving the primordial 
power spectrum with a filter function that is narrowly peaked 
at kR ^ 1 . Therefore, the abundance of UCMHs or PBHs of 
a particular mass probes the primordial power spectrum over 
a small range of scales. Meanwhile, the CMB spectral distor- 
tion produced by the dissipation of acoustic modes depends 
on the amplitude of the primordial power spectrum over a 
much wider range of scales (e.g., 50 Mpc"' ^ k < 10"* Mpc"' 
for //-distortions). Consequently, limits on the power spec- 
trum derived from CMB spectral distortions are more sensi- 
tive to the overall shape of the primordial power spectrum than 
those derived from the absence of UCMHs and PBHs. 

Possible constraints on «s and «„„ for adiabatic pertur- 
bations from future meas urements of ^ and y with PIXIE 
are discussed in detail bv IChluba et al.l ( 120121) . For a power 
spectrum with constant spectral index, PIXIE could inde- 
pendently rule out «s > 105 at ~ 5cr-level and a pure 
Harrison-Zeldovich power spectrum at ~ 2.5cr-level. This 
limit is driven mainly by the amount of small-scale power at 
wavenumbers 50Mpc"' < k < lO'^Mpc"', but a simple ex- 
trapolation from CMB-anisotropy scales down to these small 
scales is not necessarily correct. 

Here we consider more generic models for the primordial 
power spectrum that match the measurements from the CMB 
and LSS on large scales but have enhanced power on smaller 
scales. We begin in Section |2]by reviewing how the dissipa- 
tion of small-scale inhomogeneities generates CMB spectral 
distortions, and we provide A:-space window functions that fa- 
cilitate the computation of the spectral distortion generated 
by a given power spectrum. In Section |3] we compute the 
spectral distortions produced by three generic types of power 
spectra: power spectra with steps, kinks, and bumps. Most 
inflationary models that predict excess power on small scales 
generate power spectra with these types of features, and the 
constraints derived in Section |3] should be readily applicable 
to these models. We als o spe cifically discuss particle produc- 
tion during inflation (§ 13.4b . running-mass inflation (§ I3.5l l, 
and s mall -field inflation models {A<p < Mp\ during inflation; 
see § 13.61 for more details) that generate detectable gravita- 
tional waves. We summarize our analysis in Section|4] 

2. SPECTRAL DISTORTIONS CAUSED BY THE 
DISSIPATION OF ACOUSTIC MODES 

A consistent microphysical treatment of the dissipation of 
acoustic modes in t he early Universe was recently given by 
IChluba et al.l (120121) . There it was shown that temperature 
perturbations in the photon field set up by inflation lead to 
an average photon energy density (py\ =s p^[l + 6 (0^)] that 
in second order of the temperature fluctuations, = AT/T, 
is slightly larger than the energy density p^ - arT"* of 
a blackbody at average photon temperature T. The tem- 
perature anisotropics at very small scales are subsequently 
com pletely erased by shear viscosity and thermal conduc- 
tion ([Weinb erg 2008), processes that isotropize the photon- 
baryon fluid. However, the energy stored in these perturba- 
tions of the medium is not lost but merely redistributed to 
larger scales, causing a small increase of the average pho- 
ton temperature and resulting in an average spectral distor- 
tion by th e mixing of blackbody spectra with diff'erent tem - 
peratures (IZeFDovichet al.lll972t IChluba & Sunvaevll2004 . 
The effective energy release depends directly on the shape of 
the primordial power spectrum with 1/3 of the dissipated en- 
ergy sourcing y-type spectral distortions that later thermahze. 



slowly approaching a yU-type distortion. The remaining 2/3 of 
energy just causes an adiabatic increase of the average photon 
temperatu re, without crea ting any distortion. 

InChlu ba et alJ (120121) the photon Boltzmann equation for 
the average spectral distortion describing the effect of energy 
release from the dissipation of acoustic modes was derived 
and solved for primordial power spectra with constant spec- 
tral index and small runnin g using the cosmological thermal - 
ization code CosmoTherm (IChluba & Sunvaevll2012l) . It was 
shown that once the source function, (.Sac), for the primor- 
dial dissipation problem is known, a rather precise description 
of the resulting distortion can be obtained by computing the 
weighted energy release in the //- and y-era. Given (S^c) the 
required effective energy release rate caused by the damping 
of acoustic modes is determined by 



1 dgac 4f(.Sac) 



Py dz Hil+z) 



(1) 



where f = o-jN^c a; 4.4 x 10 (1 + z) sec denotes the rate 



IS 



of Thomson scattering and H ^ 2.1 x 10" (1 + z) sec 
the Hubble expansion rat^l 

Defining the visibility function for spectral distortions, 
Jbbiz) = exp(-[z/z^]5/2^, with z^ ~ 1.98 X 10^ the weighted 
total energy release in the jj- and y-era is 



Py 

Apy_ 
Py 



__ p„ ^ 

y Jo Py 



Jhh(z) dQa, 



Py dz 



dz 



dz 



dz. 



(2a) 



(2b) 



where z^,,. ~ 5 x 10^ (cf. iHu & Sill^ [19931 . At z » z^, the 
thermalization process is very efficient, so all the released en- 
ergy just increases the specific entropy of the Universe, and 
hence only raises the average temperature of the CMB with- 
out causing significant spectral distortions. However, below 
Z/j, the CMB spectrum becomes vulnerable, and energy re- 
lease leads to spectral distortions . With the simple expres- 
sions from |Sunyaev^ZeIdovic3 ([T970B), // ~ lAApy/pylfi 
and y ^ jApy/py\y; this can be used to estimate the expected 
residual distortion at high frequencies. 

The energy release depends directly on the primordial 
power spectrum, P({k), of curvature perturbations. Below we 
consider different cases for P^ik), often parametrizing it as 



P(ik) = Pf{k) + AP^ik) 



(3) 



where AP({k) describes the de viation of the power spectr um 
from the commonly used form (iKosowskv & Turrieill 19951) . 



Pf(k) = 27T^A^k-\k/ko) 



ns-l + 4"m„ln(t/*o) 



(4) 



with «,^n s dns/ d In k. In the text we often refer to Pf{k) as 
standard or background power spectrum. 

From observations with WMAP at large scales we have 
Af ^ 2.4 X 10'^ fo r pivot scale fcj_= 0.002 Mpc' 
dKomatsu et alj 1201 It iDunklev et aTj 120111 : iKeisler et al] 
120111) . Wifliout running we have «s - 0.963 + 0.014 from 
WMAP7 only, while with running the currently favored val- 
ues are «s = 1-027 + 0.051 and n,^„ = -0.034 + 0.026 

' The approximations for r and H are only valid at high redshifts during 
the radiation-dominated era. 
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jLarson et al.ll201 it iKomatsu et al.ll20T11) . More recent mea- 
surement s of the damping tai l of the C MB power spectrum 
by ACT ( lDunklevetal.ll20Tlh and SPT ( iKeisleret alj|2011h 
yielcfl ns = 1.017 + 0.036 and «„,„ = -0.024 ±0.015 and 
ns = 0.9758 + 0.0111 and «„n = -0.020±0.012,respectively. 
As Eq. ^ and our discussion below indicates, any con- 
straint derived from // or y-type distortions to leading order is 
determined by computing (i) the time average energy release 
over the redshift interval corresponding to the fi and y-era, 
and (ii) a weighted average of the total power stored over a 
particular range of scales. This means that the cosmological 
dissipation process provides a tight integral constraint on the 
power spectrum, which strongly limits possible inflaton tra- 
jectories in a very model-independent way. Furthermore, this 
constraint is not limited to just inflation scenarios but in prac- 
tice should be respected by any model invoked to create the 
primordial seeds of structures in our Universe. 

2.1. Computing the effective heating rate 

Here we are mainly interested in CMB spectral distor- 
tions caused by modes that dissipate most of their energy at 
redshifts well before the cosmological recombination epoch 
(z > 10"*), when the Universe is still radiation-dominated and 
the baryon loading R = 3pb/4-py x 673 (1 + z)"' is negligible. 
In this regime one can use the tight coupling approximation to 
compute the required source term for the photon Boltzmann 
equation (IChlubaet alJl2012[) : 

<^ac) - ^ d,k-^' f -^ k^P,(k) 2 sin^ ikr,) e-^^'I'^l, 

where ay = [1 + 4/?v/15]"- ~ 0.81 with Ry = pv/(py + Pv) ~ 
0.41 denoting the contributions of massless neutrinos to the 
energy density of relativistic species; J] - f ^ ~ 4.7 x 

10^(1 + z)"' Mpc denotes conformal time and a - {I + z)"' 
the scale factor normalized to unity today. Furthermore, 

rs a! ?7/V3 ~ 2.7 X 10^(1 + z)"' Mpc is the sound hori- 
zon; t' = ao-jNe ~ 4.5 X 10^^(1 + z)"'Mpc"' is the deriva- 
tive of the Thomson optical depth with respect to 77; ^d ~ 
4.0 X 10~^(1 +z)'^^^ Mpc~' determines the damping scale with 
(IKaiseiifTgsllZaldarriaga & Hararill 19951) 

3.9x 10^(1 +z)"^ Mpc 
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45t' 



and dimensionless sound speed c^ - Xj V3(l +R) of the 
tightly coupled photon-baryon fluid. The above expression 
for {>Sac) is based on the transfer functions for adiabatic 
perturbations, however, a similar formula can be obtained 
for isocurvature modes. S ome discussion c an be found in 
iHu & Sugivamal ( IT991 and lDent et ail ( I20T2I) . 

In the limit 7? « 1, the eflfective energy release rate for the 
photon field is therefore given by 



1 dg,, 

Py dZ 



32q' 
45 



Q-yC r dk 
Vh J 2^ 
Ckdk 
J k„ 



k^P^(k) 2 sin^ {kr,) e 



-leikl 



r^k) 2 sin^ (kr,) e 



-2k'- /kl 



(5) 



with Pik) = k^P({k)/2n^. For a given A:-mode, energy release 
happens when k ^ ku{z), where ^d is about ~ 1.9(1 + z)'^^ 

^ For both experiments we quote the constraint derived in combination 
with WMAP7 + BAO + Hq. 




Fig. 1. — Effective heating rate for the standard power spectrum, P'^Uk), 
with one sharp feature at ks- For illustration we chose (ns.ni-un) = (0.96,0). 
Furthermore, we set A? = 2.4 X 10"' in both shown cases. 



larger than the horizon scale kh ~ 2.1 x 10 ^(1 + z)Mpc ', 
implying that small-scale power is dissipated well inside the 
horizon. During the yu-era, modes with 50 Mpc ' < k < 
lO'^Mpc"' contribute most to the energy release, while y- 
distortions are mainly created by modes with A: < 50 Mpc"' . 

2.2. Energy release by a single k-mode 

We first consider the standard power spectrum with an ex- 
tremely sharp feature at some scale ks- In this case the 
modification to the power spectrum is given by APUk) - 

ItP-A^. k^^6{k - kg), where A* > determines the amplitude 
of the feature. Inserting this into Eq. (|5) yields 



1 dQiiz) 



Py 



dz 



18.8 flA^ {kslko\- sin- {ksr,) g-^A-^/*:^ ^^^ 



for the associated energy release. Notice that A:d, /"s, and a are 
all functions of redshift. This expression shows that power 
stored in a single ^-mode is released over a rather wide range 
of redshifts. The energy release peaks close to 



Zdis 



4.5 X 10' 



lO^Mpc- 



V3 



(7) 



but oscillates rapidly due to the sine part of the transfer func- 
tion. This is illustrated in Fig. [1] for ks - 8Mpc ' and 
ks = 200 Mpc"' with A^ = 2.4 x lO"** in both cases. For 

^^ = 8 Mpc"' most of the energy is released during the y-era, 
while for ks = 200 Mpc"' energy release occurs in the yU-era. 
Since the typical variation of the energy release rate is much 
longer than the oscillation period, one approximately has 



1^ ^9AaA'Uks/kufe-^'^"K 



Py dz 



(8) 



replacing sin^(x) — > 1/2, the average value over one oscilla- 
tion. The effect on the CMB spectrum can now be estimated 
by integrating the released energy over the redshifts relevant 
for the fi-eia and y-era. The y-era the integral can be per- 
formed analytically, while in the /Y-era effects related to the 
visibility of spectral distortions have to be taken into account. 
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i.e., see Eq. ^. WithyU x 1.4 Apy/py|^ andy x ^Apy/py\y the 
corresponding estimates are well approximated by 



exp - 



tis^2.2A' 



ys^OAA'lexp\- 



5400 



31.6 



- exp - 



2^ 



31.6 



(9a) 
(9b) 



where ks - kg Mpc. We mention here that the expression for 
y is only expected to be valid for kg > I Mpc"' - 5 Mpc '. 
At larger scales, baryon loading, recombination effec ts, and 
free streaming become important (IChluba et al.ll2012h . all of 
which are neglected here. These effects can be consistently 
treated using CosmoTherm, but for the purpose of this work, 
the above expression provide useful estimates for the effect on 
the CMB spectrum over a wide range of ^-values. Notice also 
that for the y-type distortions we do not apply a sharp cutoff 
in redshift, but rather limit the range in A:-space. 

2.2.1. Distortion window-function in k-space 
By replacing the amplitude A* with k^P^ik)!!?:'^ = Tik) 
and integrating over d In k, it is possible to obtain estimates 
for the values of ju and y for general primordial power spec- 
tra. Since the expressions in Eq. (|9]l are sufficiently simple, 
in many cases the integral over A\nk even becomes analytic. 
For given small-scale power spectrum we find 
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where we set ^niin - 1 Mpc"' and k = A: Mpc. These ex- 
pressions turn out to be very useful for estimates and sim- 
ple computations, as we demonstrate below. The exponen- 
tial functions act as Green's function of the cosmological dis- 
sipation problem, and the expressions can be used for gen- 
eral power spectra, as long as the effect of dissipation at 
scales k < I Mpc"' is not important. Modes in this range 
of wavenumbers are expected to affect the amplitude of the y- 
distortion, which has to be com puted using a full perturbation 
calculation ( IChlubaet al.ll2012l) . 

With these assumptions, Eq. ( fTOl i provides a weighted 
integral constraint on the small-scale power spectrum where 
fig and yg define window functions in A:-space. For a given 
detection of// this constraint has to be satisfied by any viable 
inflationary model. As we see below, COBE/FIRAS already 
placed interesting limits on several models. Furthermore, 
PIXIE will improve these limits by a large margin, strongly 
restricting possible inflaton trajectories. 

2.3. Energy release for the background power spectrum 

The total energy release and spectral distortions caused by 
the standard power spectr um, Pf{k), ac cording to Eq. (HI, 
were discussed in detail by IChluba et al.l (120 12h . with simple 
analytic approximations given for different values of ns and 
«njn. Here we are interested in cases with deviations from the 
standard shape occurring above some value of k^ni- Since the 
total ju and y-distortion are given by ju = p{k < ^ci-it) + f^{k > 
kciit) andy = yik < kcnt)+yik > k^nt), to avoid double counting 



it is therefore useful to consider the partial energy release for 
the standard background power spectrum caused by modes 
with k > kcrii- For many of our examples we s hall assume 
«mn = 0. In this case one has (IChluba et alJl2012h 

juf, ~ 5.54 X lO-'^Af exp (9.92 n^'^) , (11a) 

yf^ « 2.85 X IQ-^A^ exp (4.32 n^'^^), (lib) 

for the total fi and y-parameters. These expression were ob- 
tained using a det ailed perturbation calcula tion carried out 
with CosmoTherm (IChluba & Sunyaevl2012h . 

To compute the amount of energy release caused by modes 
with k > kcrit we start with the heating rate 
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where r(n, x) denotes the incomplete T-function. For a scale- 
invariant power spectrum we can observe the redshift scaling 
Py' dgac/ dz oc 1/(1 + z). For this reason we usually present 
the effective heating rate as (1 + z)pZ^ 'iQfc/ ^■ 

Using Eq. (|2]i one can easily compute the effective jj and 
y-parameters caused by energy release of modes k > kcdt nu- 
merically. Alternatively, with Eq. (fTOt we find 
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where k = A Mpc. These expressions work very well for 
^ciit Z 1, however, for some examples we shall use the re- 
sults obtained with a full perturbation calculation to derive 
constraints on parameters describing possible deviations for 
the standard power spectrum. 

3. SMALL-SCALE POWER SPECTRUM CONSTRAINTS 

In this section we discuss different small-scale power spec- 
tra, giving both the effective heating rates well before re- 
combination, as well as possible constraints derived from p 
and y-distortions. For COBE/FIRAS the 2cr upper limits are 
p <9x 10"^ and y < 1.5 x 10"^ dFixsen etal.ll 19961) . while 
for PIXIE one expects 2cr detection limits of yU ^ 2 x 10"^ and 
y ^ 4-x lO"** (Kogut et al. 201 1). When presenting results we 
usually assume these values, unless stated otherwise. 

In the case of COBE/FIRAS this imposes a strong upper 
bound on the amplitude of the power spectrum, while for 
PIXIE the constraints should be interpreted as 2cr-detection 
limits. Models above this limit should lead to a signal that can 
be detected at more than 2cr level, implying that they can be 
ruled out if no distortion is found. We mention, however, that 
here we do not address the more difficult challenge of using 
the detection of a CMB spectral distortion to distinguish be- 
tween different inflation scenarios. We furthermore take the 
optimistic point of view that foregrounds (e.g., due to syn- 
chrotron and free-free emission, dust and spinning dust) and 
systematics (e.g., frequency calibration, frequency-dependent 
beams) are sufficiently under control, so that the quoted de- 
tection limits of PIXIE can be truly achieved. More detailed 
forecasts including all these aspects will be required, but are 
beyond the scope of this paper. 
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3.1. Instructive upper bounds on the amplitude of the power 
spectrum at small scales 

Let us first consider the simplest ansatz for the small-scale 
power spectrum: assume that it is scale-independent with am- 
plitude !Pf (A:) - Af over some specified range of k. If we 
estimate the p and y-parameters for this case using Eq. (ITOt . 
and impose the COBE/FIRAS limits, we can determine an up- 
per bound on A^. Clearly, the constraints on A^ weaken as the 
range of scales with enhanced power narrows. 

3.1.1. Optimistic upper limit on the amplitude of the small-scale 
power spectrum from COBE/FIRAS 

Since the power spectrum at large scales is well constrained 
by CMB anisotropics and LSS observations, we first assume 
that the small-scale power spectrum has a scale-independent 
amplitude YP^ik) - 2 Ax 10 ] at all scales with wavenumbers 
k < \ Mpc"' and a different constant amplitude VP^ik) - A^] 
for A: > 1 Mpc '. In this case, one can derive an upper limit 
on Af from the p and y-limits given by COBE/FIRAS; this 
is an optimistic constraint on the small-scale power spectrum 
because we have assumed that the power spectrum is equally 
enhanced on all scales that contribute to p and y. 

Carrying out the required integrals we findQ po^ ^ 1 1 A^ 
and jop - 1.3 A^, which implies A^ < 8.4 x 10"'' and A!'' < 

1.2x10 ^ Using the TRIS bound, ;U < 6x 10"^ (IZannoni et al.l 
l2008tlGervasi et al.l2008l) . one finds A^ < 5.6x10 ^ Accord- 
ing to the weight-functions defined in Eq. ( fTOl i. the yu-limit is 
sensitive to power over 30 Mpc"' ^ k < 10"* Mpc"', while 
the y-limit is driven by the A-range 1 Mpc"' < k <, 50 Mpc"'. 
Therefore, these constraints on A^ are applicable to any power 
spectrum that has a constant amplitude over these ranges of 
scales. (Note that for y the lower cut-off is imposed by our 
assumptions, rather than for physical reasons.) 

3.1.2. Comparison to constraints from PBHs and UCMHs 

Another instructive example is motivated by the inten- 
tion to compare the power spectrum constraints derived from 
CMB spectral distortions with those obtained from PBHs and 
UCMHs. To make this comparison, we must review how the 
latter are determined from observations that limit the abun- 
dance of PBHs and UCMHs. PBHs and UCMHs form in re- 
gions where the initial density contrast exceeds some critical 
value {6plp > 0.3 for PBHs and 5plp > 10"^^ for UCMHs), 
and their masses are determined by the size of the overdense 
region that hosts them. If the perturbations are assumed 
to be Gaussian, then the probability of forming a PBH or 
UCMH with a certain mass depends only on the mass vari- 
ance within spheres that form PBHs and UCMHs with that 
mass. Therefore, an upper limit on the abundance of PBHs 
and UCMHs with a given mass implies an upper bound on 
cr^^ (/?); the density variance within a spher e of radius R eval- 
uated at horizon entry in total matter gauge (iJosan et al.ll2009l: 
[Bringmann et al. 201 1 ). These constraints on cr^^^{R) are then 
converted to constraints on the primordial curvature power 
spectrum V^ik) at wavenumber k - R\ but this conversion 
assumes a specific spectral shape for 'P({k). 

Since only the dark matter collapses to form a UCMH, 
the probability of UCMH formation depends on cr^^^ (R): 
the dark matter mass variance. In contrast, the probability 

' We confirmed these results using CosmoTherm. 



N.\ 




equivalent of PBHs 
equivalent of UCMHs - 



'"1 10 100 1000 10000 

k. in Mpc 

Fig. 2. — The bounds on A( = P^ derived from COBE/FIRAS with the same 
assumptions used to derive coiTesponding limits from PBHs and UCMHs. 



of forming a PBH depends on the total density perturbation 
at horizon entry, which is dominated by radiation. Since 
dark matter perturbations and radiation perturbations evolve 
differently as they enter the horizon, cr^^^(R) and cr^^^ (R) 
have different definitions in terms of P((k) (iJosan et al.ll2009l : 
iBri ngmann et al. 201 1). Defining x = IcR, then 



where ji (x) is a spherical Bessel function, F{x) is a filter func- 
tion, and 
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where je is the Euler-Mascheroni constant and Ci is the co- 
sine integral functi on. When e valuat ing the constraints on 
P^ik) from PBHs, IJosan et al] ( 120091) use a Gaussian filter 
function, F (x) = exp(—x^/ 2). Meanwhile. IBringmann et al] 
(2011) and iLi et"aI1 ( 120121) use the Fourier transform of a 
tophat window function, F{x) = 3.x:"''(sinx - xcosx), when 
evaluating the constraints on !Pf(A:) from UCMHs. In either 
case, if P Ak) i s nearly scale-invariant, the integrals in Eqs. 
( I14a! ) and ( I14b! ) are dominated by the contribution from a nar- 
row range of x values around x ^ 1 . To derive constraints on 
!Pf (A) from the upper bounds on cr^^^{R) and crj^^ (R) estab- 
lished by PBHs and UCMHs, it is customary to assume that 
P^ik) is locally scale invariant, i.e. that it does not vary sig- 
nificantly over the limited range of scales that contribute to 
the mass variance at a given radius. This assumption allow s 
us to take P^ik) outside the integrals in Eqs. ( |14at and ( |14bt , 
making cr^^^{R) and cr^^^ (R) proportional to P((k - 1 /R). 

Since the p- and y- distortions produced by the dissipation 
of acoustic modes receive contributions from a much wider 
range of scales than the mass variance does, there is no model- 
independent way to compare the constraints on P^ik) from 
spectral distortions to those from PBHs and UCMHs. Any 
such comparison requires one to specify the scale dependence 
of P^ik); we chose to make a comparison by applying the 
assumption of local scale invariance to the computation of 
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CMB spectral distortions. For each scale kj = l/R, we as- 
sume that P((k) is nonzero only over the range of scales that 
contribute 99% of the cr^^^{R) integral with a Gaussian filter 
(0.085 < k/ki < 2.485) and 99% of the cr^^^. (R) integral 
with a tophat filter (0.15 < k/ki < 12.35). Within these scale 
ranges, we assume that P^{k) = A^{ki) and compute the re- 
sulting spectral distortions. Since this power spectrum gives 
the same mass variance (to within 1 %) as a completely scale- 
invariant power spectrum, the constraints on the amplitude 
A^{ki) from COBE/FIRAS can now be directly compared to 
the constraints on the primordial power spectrum from PBHs 
and UCMHs; both sets of constraints make the same assump- 
tions about the local scale-invariance of the power spectrum. 

The results of our computation are summarized in Fig. |2] 
The typical limits for both the equivalent of the PBHs and 
UCMHs are Af(/t,) < 2 x lO^^ for 1 Mpc"' < kt < lO'* Mpc"'. 
At smaller scales the bound becomes less stringent because 
the thermalization process starts being very efficient. Notice 
also that for kj < 10 Mpc ' the shape of the constraint derived 
from y is affected by enforcing k > I Mpc" ' . If we omit this 
restriction the curves become practically constant at a level ^ 
10"^ for A:, < lOMpc '; however, for these cases modification 
because of recombination, baryon loading and free streaming 
should be included to obtain accurate constraints. 

For 10 Mpc"' < k < 10"* Mpc"', the upper Hmits on A((k) 
from COBE/FIRAS are more than lO"* times stronger than 
the bound obtain ed from PBHs on these scales {A^ < 0.06; 
iJosan et al1l2009l) . UCMHs can place stronger constraints on 
A({k) on these scales; if the mass of the dark matter dark 
matter particle is less than 5 TeV and it self-annihilates with 
(crv) > 3 X 10"-^cm^s"', then the fact that Fermi-LAT has 
not detected gamma-ray emission from UCMHS implies that 
Mkj) < 2 X 10"^ - 3 X lO"*" for 10 Mpc"' <k < lO^'Mpc"' 
jBringmann et al.ll201 Ih . Even if dark matter does not self- 
annihilate, UCMHs could still be used to slightly improve this 
bound; if Gaia does not detect microlensing by UCMHs, then 
A((ki) ^ 10"^ for kj ^ 3500 Mpc"' for non-annihilating dark 
matter (iLi et al.ll2012h . However, these constraints not only 
assume a particular density profile for UCMHs, but also that 
the UCMHs are not disrupted between their formation and 
today. Therefore, the bounds on A({k) from COBE/FIRAS 
are more robust. PIXIE could improve the bound on A^{k) de- 
rived from measurement of the CMB spectrum by another two 
to three orders of magnitude, potentially reaching A^ (^,) < 
few X 10^^ over scales 1 Mpc"' <k < 10"' Mpc"'. 

3.2. Constraints on steps in the power spectrum 

Next, we consider a step in the power spectrum at some 
scale k^, where the amplitude changes from A^ for ^ < A:s to 
Af + A^ for k> ks- Such a step in the primordial power spec- 
trum could be produced by multi-stage inflation models or in- 
flaton potentials t hat change slope when the inflaton reaches 
a certain value (ISilk & Turned 119871: ISal opeket all 119891 
Starobins kiill992l;lPolarski & S tarobinskv 1 992llIvanovet all 
1994: .Adams et al.lll997t IStarobinskvMl998lK We will as- 
sume that the spectral index is the same on both sides of 
the step, and we note that the step has to fulfill the condi- 
tion AJ. > -Af , as otherwise unphysical negative power is 
present in the power spectrum. These power spectra can be 
parametrized by Eq. d^i with ^P^ik) = A^P^'C^VAf at A: > Ats 
and Af ^ - otherwise. If for simplicity we assume «„„ - 
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Fig. 3. — Effective heating rate for the standard power spectrum, P^}(k), with 
a step at ij. For illustration we chose («s. "run) = (0.96, 0). Furthermore, we 



set A^ = 3.5 X 10 'in both shown cases. 
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Fig. 4. — Limits on the total amplitude of one step in the small-scale power 
spectrum occurring at is. The energy release of modes with k < k^ was 
included, however, this only matters for the case of PIXIE. For illustration, 
the power law index was also varied. 



for the background power spectrum, P^', we find the efifective 
heating rate by modes with k>ksis simply given by Eq. ( fTSl l 
with A^ replaced by Aj^+A\. Similarly, the ^ and y-parameters 
caused by the change in power can be estimated using the ex- 
pressions Eq. ( fT3T l. 

In Fig. [3] we illustrate the time-dependence of the effective 
energy release for k^ - 10 Mpc"' and k^ - 200 Mpc"' with 
step amplitude A^ - 3.5 x 10"^ in both cases. The total 
amplitude of the power spectrum after the step therefore is 
A'°' - A( + A!^^ 5.9 X 10"''. In contrast to the single-mode 
case, we see that the energy release no longer exhibits any os- 
cillatory behaviour, since oscillations of several neighbouring 
modes cancel each other, leading to smooth energy release. 

In Fig. ID we present constraints on the total amplitude 
Af + Al of the power spectrum after the step. We show 
both limits obtained from COBE/FIRAS and possible future 
bounds from PIXIE. In the COBE/FIRAS case the energy re- 
lease from the background spectrum, Pf{k), can be neglected, 

as it only results in // ^ y ^ 10"^ for reasonable values 
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of «s and «„„ ( IChlubaet al.ll2012h . The COBE limits ob- 
tained from the /^-distortion are most stringent in the range 
ks < lO-' Mpc"'; the Hmits get weaker at smaller scales be- 
cause the thermalization process starts being very efficient. 
Also the /i-limit is stronger than the y-limit because the fi dis- 
tortion receives contributions from a slightly larger logarith- 
mic range of wavenumbers (1 Mpc"' i k < 50Mpc"' fory as 
opposed to 50Mpc"' < A: < lO'^Mpc"' forju). 

For PIXIE we only present the possible constraint derived 
by measurement of fj.. Obtaining a limit from y is expected 
to be much harder, because at low redshifts many other as- 
trophys ical processes ( e.g., energy release because of super- 
novae JOhet alj 12003b: shocks dur ing large scale structure 
formation ( Sunvaev & Zeldovichll9 72; Cen & Ostriker 1999; 
Miniati et al.ll2000l) ; unresolved SZ clusters ( iMarkevitch et akl 
19911) : the thermal SZ effect and second ord er Doppler ef- 
fect from reionization (iMcOuinn et al.l 120051) ) can cause an 
average y-distortion that is expected to be orders of magni- 
tude larger in amplitude. For PIXIE the energy release caused 
by the background power spectrum no longer can be ignored, 
since fi(k < k^,) ^ fi{k > k^). We therefore present the limit 
on just the amount of dissipation at scales k > k^, after sub- 
tracting fi{k < k^) as given by Eq. (fT3] l. For the shown ex- 
amples with ns - 1, «s = 1.05 and «s = l.I our results 
imply that respectively for k^ > 100 Mpc"', 250 Mpc"' and 
10 Mpc"' a more than 2cr-detection of yu-distortions is ex- 
pected even in the case A^ + A^. = for k > k^. This is simply 
because fi{k < k^) itself already exceeds the 2cr-detection limit 
of PIXIE, i.e., //lin, ^ 2 X 10^1 Also, at about k ^ 20Mpc-' 
the curves become flat, indicating the point at which practi- 
cally no additional //-distortion is produced by modes with 
smaller wavenumber In this regime only the amplitude of 
the y-distortion is expected to change; however, unless a large 
distortion (y > 10"^) is created, this signal will be hard to sep- 
arate, as mentioned above. Nevertheless, simultaneous detec- 
tion of (large) y and jj could constrain the scale at which the 
step occurred. 

Here we only considered one step, but it is easy to general- 
ize the discussion to multiple steps. The bound will strongly 
depend on the distribution of k^^ and Al for which physical mo- 
tivation should be provided, suggesting a case-by-case study 
is more useful. If, for example, all steps increase the total 
power at small scales, then the derived limits are expected 
to become stronger However, for models with oscillations 
around a standard scale-invariant small-scale power spectrum 
(see, e.g., lAchucarro et al. 2011: .Cespedes et al..,2012,) . the 
net effect should average out. 



3.3. Constraints on a bend in the power spectrum 

As a second example we consider a kink or bend in the 
power spectrum at some scale ^b with the slope of the power 
spectrum changing from «s to «g, while P{k) remains con- 
tinuous. iJov et al.l (|2008) showed that such changes in the 
spectral index result from discontinuities in the second deriva- 
tive of the inflaton potential, and this model for P{k) may 
be used to approximate the power spectra produced by sev- 
eral other m odels tha t genera t e large perturbations on small 
scales (e.g.. IStewarj Il997at lBen^Davan& Brusteinl 12010 : 
Gong & Sasaki"201 1; 'Lyth"201 la"; 'Bugaev &Kliinai|l2011a|; 
Shafi & W ickman,i201 k iHotchkiss et al., i2012l) . The associ- 
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Fig. 5. — Effective heating rate for tlie standard power spectrum with a kink 
at^b. For illustration we chose (ns,«run) = (0.96,0). Furthermore, we set 
«j = 1,2 in both examples. 



ated power spectrum can be parametrized as 

* |27r2r^V^,'(fcb)(/t//tb)"s-' at/t>fcb- 



(16) 



with !P^'(A:) = k^ Pf(k)/27T-. Assuming «„,„ = 0, the to- 
tal power released by modes with k > kh is again given by 
Eq. ([HI with A( replaced b}0 A^ (kb/ko)"''-"s. Similarly, the 
// and y-parameters caused by modes with k > k\, can be esti- 
mated usingthe expressions Eq. (flJl l. 

In Fig. \5\ we illustrate the effective heating rate for 
(«s,nnjn) - (0.96,0) and n^ = 1.2. One can clearly see a 
flaring of the energy release that starts close to Zdiss(^b) ac- 
cording to Eq. dTji. The redshift dependence suggests that the 
effective y-parameter caused by a bend in the power spectrum 
is typically smaller than the //-parameter Indeed we find that 
for the COBE/FIRAS limits the constraints derived from y 
are much weaker than those from //, so we neglect them for 
the discussion below. Once again, concerns regarding confu- 
sion with low redshift y-distortions are the limiting factor for 
constraints derived from PIXIE's measurement of y, although 
simultaneous detection of both // and y from the cosmological 
dissipation process would provide a deeper understanding of 
the shape of the small-scale power spectrum and the position 
of a possible kink. 

We can derive constraints on the value of «g, and these 
constraints are shown in Fig. |6] The COBE/FIRAS limit 
on fi already rules out changes in the power law index by 
Ans = Hg - «s ^ 1 at ^b - 1 Mpc"' at 2cr-level. With 
PIXIE this measurement will be strongly improved. For ex- 
ample, if the background spectrum has ns - 0.96 then even 
at k\, ^ 100 Mpc"' the slope cannot change by more than 
A«s - 0.1 without leading to an observable //-distortion. 
Also, like in the step case, for «s ^ 1 a bend at large values 
of ^b should lead to an observable signal even if the power 
spectrum cuts off abruptly (i.e., «g has a large negative value). 
This should allow placing very tight constraints on inflation- 
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Fig. 7. — Effective heating rate for the standard power spectrum with a fea- 
ture at kj caused by particle production during inflation. For illustration we 
chose («s,'!nin) = (0.96,0) for the background spectrum. Furthermore, we 
set A^' = 3.5 X 10 in both shown cases. 



ary models with flaring power spectra at small scales. 

We note that PIXIE cannot constrain hybri d inflation mod- 
els th a t use a 'wate r fall' field to end inflation (iGong & Sasakil 
120111; lLvflj|2011at iBugaev & Khmail2011al) because these 
models predict k^, '» 10'^ Mpc"', corresponding to scales that 
left the horizon during the last few e-folds of inflation. Even 
PIXIE will be insensitive to energy release from those scales. 

3.4. Constraints on particle production during inflation 

iBamabv et all (l2009h and iBarnabvl ( 120101) showed that 
bursts of particle production during inflation produce local- 
ized bumps in the primordial power spectrum. Specifically, 
if the inflaton is coupled to another scalar field x via an 
interaction given by Xim = -\g]{<P " <Pif'X^^ then the;^' par- 
ticles are temporarily massless when (p - (p,. At this time, 
X particles are created by quantum effects, and these parti- 
cles quickly become massive as (p moves away from 0,. The 
massive x particles then rescatter off the field, generating 
perturbations in that freeze once their wavelength exceeds 



the Hubble distance. The massive x particles are rapidly di- 
luted by the inflationary expansion, so only a limited range of 
s cales receive ex tra perturbations. 

iBamaby & Hua ng (2009) provided a simple parametriza- 
tion for the the resulting bump in the primordial power spec- 
trum: 



ApP''(fc) 



2;r-AP'' 



3/t? 



3/2 



^^P|--72 



1^ 

'2k 



(17) 



The amplitude of the feature, A^'', is simply related to the 
value of coupling constant gi\ A^'' ^ 1.01 x 10 "g'. 



6„15/4 



The 



derivation of this feature in t he primordial power s pectrum is 
only valid for 10"^ ^ gf ^ 1 ( iBamabv et al.ll2009b . so we are 

primarily interested in A^' values between 10 '^ and 10"^. In 
contrast, there are no restrictions on the location of the bump; 
kj is determined by the number of e-foldings between the mo- 
ment when (p - (pi and the end of inflation. There may also 
be other fields with the same coupling to the inflation, each 
with their own values for g, and (p,. In this case, the power 
spectrum will contain multiple bumps, and one should sum 
the contributions from each episode of particle production. 
Inserting Eq.fTTlinto Eq. (|5j, we find 



1 dgp,, 



Py 



dz 



9AaA 



(ki/kuf 



( 



Vfo^ [ 



1 + 



'^(ki/koy^] 



5/2 



(18) 



for the effective heating rate of one feature at high redshifts 
(see Fig.|7]for illustration). The fi and y-parameter caused by 
one feature are roughly given by 



;Up - 2.4 A^A 



e(f{l+2^^)ETfc(^i)- 



2^i yp 



V^ 0.4 A^-'^ 



11-3/2 



(19a) 
(19b) 



3'p^0.4AP'' [l+(fc,/40)- 

where ki = ki Mpc, ^; - ki/1.35 x 10"*, and Erfc(x) is the com- 
plementary error function. We again made use of Eq. (ITOi 
to give these simple expressions, but we also confirmed the 
validity of these expressions by numerically evaluating the 
nested integrals of the power spectrum and the heating rate. 

In Fig. [8] we present the derived limits on the amplitude for 
one episode of particle production. The limits derived from 
COBE/FIRAS are weaker than the bound A^' < 10 *> that is 
required to make the underlying calculation self-consistent. 
However, these bounds are still interesting, as they can be 
also interpreted as rather tight constraints on any other infla- 
tion models with bump-like features in the small-scale power 
spectrum that have a typical total width of Ak/k ~ 2/ V3. 

The bound becomes much tighter for PIXIE, basically lim- 
iting A^'' < 2x10-'' in the range lO'Mpc-' <k< lO^^Mpc"'. 

Also, in the range 1 Mpc ' < k < 100 Mpc ', the hmit 
derived from ju alone is still very interesting, although it is 
weaker since fewer modes related to the bump are able to re- 
lease energy in the yu-era. Similarly, the bound becomes less 
stringent for k > lO"* Mpc"' because thermalization becomes 
efficient. Features in the small-scale power spectrum intro- 
duced by particle production are rather broad, with a signifi- 
cant tail of energy release towards lower redshifts, where the 
visibility for spectral distortions increases. Therefore, even 
for larger ki an observable /i-distortions is created, soften- 
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Fig. 8. — Limits on inflation models with particle production. Here we 
only considered one episode of particle production. Tiie dasiied line indi- 
cates tire bound a"!' < 10"'', required to ensure physically reasonable models 



jBamabv & Huang|2009f) . We assumed a background power spectrum with 
(ns,nrun) = (0.96,0). 
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Fig. 9. — Effective heating rate for inflatio n m odels with running-mass. The 
power spectrum is pai'ametrized as in Eq. I I20I . For reference we also show 
the standard power spectrum for («s, "run) = (0.96, 0). 



ing the thermalization cut-ofF. This implies that PIXIE could 
even constrain episodes of particle production with k\ up to 
~ 5 X 10"^ Mpc"' in an interesting way, complementing the 
limits obtained from CMB an isotropies and LSS at larger 
scales (iBamabv & Huangll2009l) . 

We also mention that for more than one episode of parti- 
cle production, the constraints should become tighter How- 
ever, in this case again physical motivation should be given 
and models ought to be discussed on a case-by-case basis. 
Furthermore, the constraint on bumps in the power spectrum 
also depends on the background power spectrum at small 
scales. For models with ns - 1 direct detection of the yi- 
distortion with PIXIE could be possible even without extra 
bumps, so that any excess power added by particle production 
features should further enhance the yt/-distortion above the 2cr 
detection threshold. This indicates that the interpretation of 
the constraint depends significantly on assumptions about the 
background model extrapolated from CMB and LSS scales all 
the way Vok^ lO'* Mpc"'. 



3.5. Constraining running-mass inflation models 

For our next explicit example, we consider running-mass 
inflation, which is a single-field supersymmetric inflation 
model that can generate enhanced power on smal l scales 
(lStewarlll997allbl:lCovi &. Lvflill999tlCov"i et al.ll999h . These 
models assume that the dominant loop correction to the infla- 
ton potential is oc 0^1n[0/g], where Q is the renormalization 
scale. In this case, the renormalization-group-improved in- 
flaton potential is y((^) = Vq + m^(\n^)^- jl, implying that 
the mass of the inflaton effectively changes during inflation. 
Consequently, running-mass models offer a solution to the 77- 
problemofsupergravity inflation; in supergravity, scalar fields 
usually have masses that are too large to drive inflation, but in 
running-mass inflation, \m^\ can be large when (p is equal to the 
reduced Planck mass Mp\ while being small enough to permit 
inflation at smaller values of (p. 

To ensure that the inflaton potential is sufficiently flat, in- 
flation must occur near an extremum of the inflaton potential. 
Therefore, we can approximate the inflaton mass as 



2 c Vo 
pi 



1 



21n — 

0. 



where (pt is the value of at the extremum of y(0) and 



-(M^./yo)dm^/dln0 evaluated at 



The 



parameter may be positive or negative, but explicit formula- 
tions of running-mass i nflation in the co ntext of supersymme- 
try have < c < 1 (e.g. lCovi et al.l2004l) . Given this linear ap- 
proximation for n-p-jln S), the primordial power spectr um can 
be parametrized as (cf. lCovi & L"vflilll999l:ICovi et al.ll2004 



nfc) = !P(^o)exp 



2s 



ii(e'-AW_i)_2cAA^(;t) 



(20) 



In this expression AN(k) = ln{k/ko), and s - c \n{4>t 1 4>(i)^ 
where </>() is the value of the inflaton when the mode with 
wavenumber k^ exited the Hubble horizon, A^o e-folds before 
the end of inflation. Note that s can be positive or negative, 
depending on the sign of c and the direction (p rolls during in- 
flation. It follows that s - cln((^»/(/ie)exp[-cA^o], where 0e is 
the inflaton value at the end of inflation. Therefore, values of 
s with 1^1 < |c|e"^''° imply that t^e - 0*, which requires certain 
fine-tuning. 

One can also directly relate s and c to the usual spectral 
index and running dCovi e t al. 2004): 



ns ^ 1 + 2 [s{l - 1.06c) - c] + - (.? - cY , (21a) 



2sc. 



(21b) 



Observational parameter limits on s and c derived usi ng pre- 
vious WMAP (Bennett etal. 2003), SDSS (TegmarlmaD 
120041) and Ly-o' forest (ISeljak et al.l 2005) measurements are 
shown as closed contours in Fig. [TOl (Covi et al...2004.) . lim- 
iting the range of allowed models to cases with s =^ c and 
\s,c\ < 0.13. According to Eq. (ISTT i. the allowed models typ- 
ically have positive running with -0.001 < n^n S 0.02 and 
0.91 < «s ^ 106. This indicates that updated constraints 
from the latest CMB and LSS measurements that favor neg- 
ative running might already further narrow down the allowed 
parameter space in comparison to ICovi et all (120041) . How- 
ever, the two running-mass models that match the current 
best-fit WMAP7 model without running, (ns, «i™) = (0.96, 0) 
at fe = feo are still viable (cf. Fig.fTOt. 
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Fig. 10. — Constraints on ranning-mass inflation models obtain b y measure- 
ments of//. The closed contours show the constraints obtained by jCovi et aH 
I2004D . PIXIE might help naiTowing down the allowed pai'ameter space to a 
small region around .v =: c - 0.033 with -0.03 < c < 0.06. Furthermore. 
PIXIE with twice the sensitivity could in principle rale out running-mass in- 
flation models in case no distortion is found. For comparison we also show 
the location of the two models that are consistent with (ns.nrun) = (0.96,0) 
at the pivot scale io- 

Using Eq. (|20] i we can easily compute the effective heating 
rate and distortion parameters for different values of s and c. 
In Fig.|9]we illustrate this for two runnin g-mass models tha t 
are in agreement with the constraints of ICovi et al.l (120041) . 
The departure from the standard background spectrum occurs 
very gradually, so that in comparison with the standard power 
spectrum, both enhanced /i and y-distortions are expected in 
nearly all cases. 

To compute the limits on the parameters s and c shown in 
Fig. [To] we use the expressions given by Eq. ( fTOl i. We also 
confirmed that a more precise computation practically gives 
the same result. The COBE/FIRAS limit on /i does rule out a 
significant part of the theoretically allowed param eter space, 
however, in comparison to the constraints given bv lCovi et al.l 
(1^04) no improvement is achieved. On the other hand, PIXIE 
might rule out many running-mass inflation models at the 2cr 
level, narrowing the possible parameter space down to a slim 
region around s ^ c - 0.033 and -0.03 < c < 0.06. Further- 
more, PIXIE with twice the sensitivity could rule out running- 
mass inflation models if no distortion is detected. We also 
computed the limit from measurement of the y-parameter, but 
found that the constraint was always much weaker 

3.6. Spectral distortions for small-field models generating 
detectable gravitational waves 

Another interesting class of inflation models that predicts 
an enhancement of the power spectrum at small scales is a 
class of small-field models th at generate a detectable grav- 
itational wave (GW) s ignal ('Ben-Da van & BrusteinI 120101 : 
IShafi & Wick"rnanll20rit [Hotchkiss et alj|2012t). wifliin this 
class, models with larger GW signal are also expected to pro- 
duce excess power at small scales. Therefore, the detection of 
both a GW signal and a CMB spectral distortion at the pre- 
dicted level would provide strong evidence for these models. 

The exa ct power spectrum has to be e valuated numerically 
(as done in lBen-Davan & Brusteinl201CII) and exhibits a richer 
phenomenology than the standard nearly scale-invariant case. 



even when «run is included. The ability to produce a detectable 
GW signal with limited field excursion, tS.(p < Mpi, is based 
on the fact that the slow-roll parameter e is non-monotonic 
from the time the CMB scales left the horizon until inflation 
ends. More explicitly, e is initially large enough to produce a 
detectable GW signal, then it decreases to a smaller value for 
most of the duration of inflation. At a later time, e increases 
to unity and inflation ends. The scales of 10 "* Mpc ' < k < 
lO"* Mpc ' exit the horizon while e is decreasing; since !Pf oc 
y/e, small-scale power is enhanced. The power spectrum also 
exhibits a maximum at a smaller scale that exits the horizon 
just before e begins to increase. 

Since here we are only interested in the possible CMB dis- 
tortions we do not repeat the exact numerical calculation, 
but use the 'improv ed power spectrum' form (Eq. (78) of 
iLyth & Riottolll999h . which captures all essential features of 
the model. A more complete likelihood analysis is deferred 
to later work. We c alculated the n and y dis t ortion s for some 
models discussed in lBen-Davan & BrusteinI (120101) . however, 
we added specific model^ with small tensor to scalar ratio 
r ^ 10 "'. The GW signal of these scenarios could potentially 
be detected with PIXIE. 

As mentioned above, we find that the level of /i and y 
distortions is linked to the GW signal. This is because a 
larger GW signal requires a larger change in e between CMB 
scales and smaller scales. We therefore expect a bigger 
enhancement of small-scale perturbations. For one of the 
models with relatively large GW signal, r = 0.08, we obtain 
fi ^ 8.5 X 10"^ and y ^ 1.3 x 10"^. At this level, we expect 
a definite detection by PIXIE. However, the distortions 
decrease with r, so that for r ^ 10"^ the spectral distortions 
reduce to the level similar to the standard po wer spectrum 
without running, i.e., fu ^ 10"^ and y ^ 10"^ dChluba et al.l 
120121) . Therefore, we have an interesting cross-check be- 
tween GW detection and CMB distortions. In the optimistic 
case, if the model of inflation realized in nature is of this type, 
then we expect PIXIE to detect both tensor perturbations and 
CMB distortions. 



4. DISCUSSION AND CONCLUSION 

We considered constraints on the small-scale power spec- 
trum derived from present and future measurements of the 
spectral distortions in the CMB. We introduced A:-space win- 
dow functions for ju- and y-distortions that account for the ef- 
fect of thermalization and dissipation physics and facilitate 
computing the effective fi- and y-parameters directly from the 
primordial power spectrum at ^ > 1 Mpc"' (see Eq. ( fTOt in 
§ 12.21 ). This defines an integral constraint that places tight 
limits on the shape of the small-scale power spectrum and, by 
extension, constrains possible inflation scenarios or any other 
model of the early Universe that is invoked to create the pri- 
mordial seeds of the structures we see today. 

We discussed different generic cases for the small-scale 
power spectrum, demonstrating how this integral constraint 
can be translated into limits on power spectrum parameters. In 
particular, we derived limits from the COBE/FIRAS bounds 
on fi and y, showing that for 1 Mpc ' < k < 10"* Mpc ' 
these upper limits on the amplitude of the power spectrum 
are roughly 10^ times s tronge r than those derived from PBHs 
at similar scales (see § 13.1.21 ). Limits obtained with UCMHs 

' Following the definitions of lBen-Davan & BrasteinI 1201(1) . these corre- 
spond to: ffo = ;;o = and i^end = 1 ; co = 0, ;;o = -0.02 and i^end = ' ■ 
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supersede the COBE/FIRAS limits by more than an order of 
magnitude, but these limits depend on the properties of the 
dark matter particle. In contrast, the constraints from CMB 
distortions can be obtained in a very model-independent way, 
relying on well-understood physics related to the thermal- 
ization and dissipation of acoustic modes. We also showed 
that PIXIE will improve the bounds derived from ji and y 
by about three orders of magnitude. PIXIE could therefore 
open a new window to the early Universe, extending the lever 
arm from CMB-anisotropy and LSS scales all the way to 
k^ lO^Mpc-' - lO'^Mpc-'. 

As explicit examples, we studied the constrai nts on inflation 
models with episode s of particle production (§ I3.4| | and run- 
ning inflaton mass (§ 13. 5b . We demonstrated that PIXIE could 
complement the upper limits on particle production derived at 
CMB and LSS scales, extending them from ^ < 1 Mpc ' up 
iok ^ 10"* Mpc ' (see Fig. ID . We also showed that PIXIE 
might have the opportunity to rule out running-mass inflation 
models if no spectral distortion at the level of // ~ 2 x 10"'^ is 
found (see Fig.[TO]and § l3.5l l. Similarly, for other models with 
flaring small-scale power spectrum (e.g. iBarnabv et alll201 lb 
our computations indicate that strong bounds could be placed 
on the viable parameter space. As argued in § 13.61 small-field 
inflation models {A(f> < Mp\ during inflation) with significant 
GW signal should simultaneously produce large CMB spec- 
tral distortions, an intriguing connection that could potentially 
be established by PIXIE. 

The possible limits derived from future PIXIE measure- 
ments of y suffer from confusion with y-type distortions cre- 
ated at low redshifts because significant differences between 
these two signal are not expected. However, one effec t migh t 
help in this respect: as shown by IChluba & SunvaevI (l2009h 
energy release before recombination causes uncompensated 
cycles of atomic transitions in helium and hydrogen. This 
leads to extra emission and absor ption features in the cos- 
mological recombination radiation (IChluba & Sunyaevll2006l: 
ISunvaev & Chluball2009l) . At both high and very low fre- 
quencies, the associated effect is larger or comparable to the 
y-distortion itself, so that a delicate interplay between the 
thermalization process and the recombination radiation is ex- 
pected. This effect might distinguish y-type distortions im- 
printed before the end of recombination from those coming 
from lo w redshifts, if precise spec tral measurements are per- 
formed (IChluba & Sunvaevll2009l) . In particular, distortions 
at high frequencies, created by the Lyman-series of hydrogen 
and doubly ionized helium, as well as the n ' P i - 1 ' So series of 
neutral helium might be very interesting in this respect. How- 
ever, to give a definite answer a more detailed computation is 
required, simultaneously including the effect of atomic transi- 
tions in the thermalization calculation. 

The derived bounds are obtained under the assumption that 
the only process causing energy release at high redshifts is the 
dissipation of acoustic modes. However, other forms of en- 
ergy release are pos sible. For example, as r e cently shown by 
IChluba & SunvaevI (12012b and iKhatri eraLJ (120111) . the adia- 
batic cooling of ordinary matter inevitably leads to small neg- 
ative ji- andy-type distortions with amplitude /i ^ -2.4 x 10"'' 
and y ^ -4.3 x 10 "I This process is based on well un- 
derstood physics, and hence the associated distortion can be 
predicted with high precision. Nevertheless, many other pos- 



sible sources of energy relea se exist including amiihilating or 
flfecay/ w g particles (e.g., seelBuriganaet al.lll99lt iHu & Silfl 



1993allb :lMcDonald et al.ll2001h : evaporatins black holes (see 



Carret 



y_ 20101 and references therein); superconduc ting 



strings (lOstriker & Thompsonl '1987: Tashiro et al. 20121); or 
dissipation of magnetic f ields ( Jedamzik et al. 2000) . (See 
IChluba & SunvaevI (120 12b for detailed computations of the as- 
sociated distortions with CosmoTherm.) 

All these processes come with significant uncertainties, so 
it is unclear at which level distortions can be expected. There- 
fore, the constraints obtained here should be considered as 
the most conservative upper limits, since any additional en- 
ergy release not caused by the dissipation of acoustic modes 
will only tighten the bounds on the primordial power spec- 
trum. An interpretation of a detection of CMB spectral distor- 
tion therefore requires more careful consideration of the dif- 
ferences (e.g., the mixture between// andy; the detailed shape 
of the distortion at low frequencies) in the distortions for each 
case, which in principle can be accurately computed using 
CosmoTherm. In addition, differences in the spatial distribu- 
tion, although expe cted to be tiny for primordial distortions 
( IChluba et al.ll2012h . might help distinguish different scenar- 
ios in the future. Correlations of the distortion with CMB 
anisotropics coul d furthermore reveal non-G aussianity of the 
power spectrum (iPaier & Zaldarriagall2012b . 

We close by mentioning that even if the parameters 
describing the power spectrum at CMB-anisotropy and 
LSS scales fully determine the small-scale power spec- 
trum, a measurement of /i and y could be interpreted as 
an independent confirmation of these values. Moreover, 
it is important to note that the determination of «s and 
«i-un with CMB measureme nts is subject to unc ertainties 
in recombination dynamics (IShaw & Chluball201 lb . While 
standard recombination physics seems to be under control 
(e.g., see Dubrovich & Grachev 2005; Kholupenko et al] 
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possible surprises due to neglected standard or non-standard 
processes could still await us. Directly constraining the 
recombination dynamics fr om CMB anisotropy measure- 
ments itself is challenging (iFarhang et al.l 12011b . so some 
theoretical uncertainty in the values of «s and «nan is un- 
avoidable. Therefore, a detection of yu-distortions at the level 
extrapolated from CMB and LSS scale constraints would be 
very reassuring, further demonstrating the great potential of 
this new window to the early Universe. 
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